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Abstract 

We investigated numerically an Ising model coupled to two- 
dimensional Euclidean gravity with spherical topology, using Regge 
calculus with the dl/l path-integral measure to discretize the gravita- 
tional interaction. Previous studies of this system with toroidal topol- 
ogy have shown that the critical behavior of the Ising model remains in 
the flat-space Onsager universality class, contrary to the predictions of 
conformal field theory and matrix models. Implementing the spherical 
topology as triangulated surfaces of three-dimensional cubes, we find 
again strong evidence that the critical exponents of the Ising transition 
are consistent with the Onsager values, and that KPZ exponents are 
definitely excluded. 
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1 Introduction 



The study of two-dimensional (2D) models has often proved to be an im- 
portant first step for developing fundamental ideas concerning higher di- 
mensional physics. Such has happened for 2D Euclidean quantum gravity 
where we have by now independent analytical results from conformal field 
theory and matrix models [1]. For pure gravity this is a formula for the 
string susceptibility exponent 7 str , and for matter fields of central charge 
c < 1 coupled to gravity it has been found that the nature of the phase 
transition changes drastically. In the simplest case of the Ising spin system 
(c = 1/2), the critical exponents change from the flat-space Onsager values 
a = 0, = 0.125, 7 = 1.75, and v = 1 to the KPZ [2] values a = -1, 
(3 = 0.5, 7 = 2, and Du = 3, where D is the internal fractal dimension of the 
manifold. While these predictions could be confirmed numerically within the 
dynamically triangulated random surface approach, it was found [3, 4] that 
the Ising model on lattices of toroidal topology stays in the flat-space On- 
sager universality class when coupled to gravity via Regge calculus [5]. This 
came as a big surprise because on the basis of universality one had expected 
that Regge calculus, as the best understood method to discretize classical 
gravity, would give results compatible to KPZ scaling. The Regge method 
seems to depend sensitively on the path-integral measure [6], and therefore 
the prevailing opinion is that the measures used so far are inappropriate [7]. 
This issue is, however, by no means settled and still requires further investi- 
gations. It is also not clear if the commonly chosen spin coupling to gravity 
is the correct one. 

Furthermore, in pure gravity the string susceptibility exponent 7 str de- 
pends on the global topology, and one might wonder if the spin phase tran- 
sition also shows a dependence on the global topology. In particular the 
previously used torus topology plays a special role, since 7 str takes the classi- 
cally expected value of two, suggesting the absence (or triviality) of quantum 
effects in this topology. In order to exclude the possibility of a fortuitous co- 
incidence with flat-space critical exponents in the torus topology, we have 
therefore performed further Monte Carlo simulations of the Ising system us- 
ing the lattice topology of the sphere. 
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2 The model and simulation 



The model and our notation are explained in detail in Ref. [4], where we 
used a regular triangulation of a torus. We recapitulate here only the essen- 
tial ingredients. Because in two dimension the Einstein-Hilbert action is a 
topological invariant, we need no curvature term. In Ref. [4] the Ising system 
was studied with various curvature squared terms, which showed no influ- 
ence on the critical behavior of the system. Therefore in the present study 
we chose to have no R 2 term at all. This leaves only the integral over the 
metric g to be discretized, which gives the total area A. We used 

i 

where the A{ are taken as the barycentric areas associated with the vertices 
i. We simulated the partition function 

Z = J2 [ exp (-X Ai - KE(l, s)) A, - A), (2) 

where 

E(l,s) = \ £ (3) 

edges (ij) 

is the energy of Ising spins, s 8 - = ±1, which are located at the vertices i of 
the lattice, and Aij is the barycentric area associated with an edge (ij). The 
energy is the discretized analogue of the continuum action for a scalar field 
<f, J cPx^Jg(x)g tJ,u (x)d fl <p d v ip. The delta function in (2) ensures that the total 
area is kept fixed at a given value A. As path-integral measure T>fi(l) we used 
again the simple scale invariant measure T>/i(l) = (jl^j) dUj/U^j F({lij}), 
where F({lij}) is a function which ensures that changes in the link lengths 
do not violate the triangle inequalities. 

As lattice topology we used the triangulated surface of a three dimensional 
cube of edge length d. This provides us with an almost regular triangulation 
of the sphere where six vertices have coordination number four, and all others 
have coordination number six. In terms of the linear length d of the cube 
the number of vertices is Ao = 6(o? — l) 2 + 2. For further reference the 
number of links and triangles in terms of Ao are given by A^i = 3Aq — 6, 
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and N 2 = 2N — 4, respectively. We studied ten system sizes ranging from 
d = 10(iVo = 488) up to d = 55(N = 17498). The area was kept fixed 
to its initial value A = N 2 /2. The gravitational action was simulated using 
the standard single-hit Metropolis update. For the spin update we used 
the single-cluster (Wolff) algorithm [8] which prevents the critical slowing 
down near the phase transition. Between measurements we performed n = 2 
Monte Carlo steps consisting of two lattice sweeps to update the link lengths 
lij followed by one single-cluster flip to update (a fraction of) the spins 
As simulation point we chose a value of K = 1.025, which is close to the 
critical coupling found in Ref. [4]. 

For each run we recorded the time series of the energy density e = E/N 
and the magnetization density m = J2i Ai-Si/N . After an initial equilibration 
time, we performed for each lattice size about 50 000 measurements. From an 
analysis of the time series we obtained integrated autocorrelation times for 
the energy and the magnetization of about 1 — 7 (in units of measurements) 
for all lattice sizes. The statistical errors were obtained by the standard 
Jack-knife method using 20 blocks. 

From the time series we computed the Binder parameter, 



apart from confluent corrections explaining small systematic deviations. This 
allows an almost unbiased estimate of the critical coupling K C} the critical 
correlation length exponent z/, and the renormalized charge U*. We further 
analyzed the specific heat, 




(4) 



for each lattice of size L = y/N , which serves as our linear scaling parameter. 
The curves Ul(K) for different L cross around (K C ,U*) with slopes oc L 1 ^, 




(5) 



the (finite lattice) susceptibility, 




(6) 



the susceptibility in the disordered phase, 




(7) 
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and studied the (finite lattice) magnetization at its point of inflection, (|m|) |; n f . 
The inflection point can be obtained from the maximum of d(\m\) / dK } which 
can be calculated from the time series as 

*^ = (E)(\m\)-(E\m\). (8) 

We further analyzed the logarithmic derivatives dln(\m\) / dK and o?ln(m 2 } / dK . 



3 Results 

By applying reweighting techniques we first determined the maxima of C, 
X, d(\m\) / dK , d\n(\m\) / dK , and dln(m 2 ) / dK . The location of the max- 
ima provided us with five sequences of pseudo-transition points K max (L) for 
which the scaling variable x = (K ma , x (L) — K c )L l l v should be constant. Using 
this information we then have several possibilities to extract the critical ex- 
ponent v from (linear) least square fits of the finite-size scaling (FSS) Ansatz 
dllj^ldK = L l l u 'fo(x) or dln(\m\' p ) / dK = L 1 ^ 1J f p (x) to the data at the various 
K ma , x (L). The resulting exponents 1/u from fits with an acceptable goodness- 
of-fit parameter Q are collected in Table 1. They give mostly results which 
are smaller than the Onsager value v = 1, but are still compatible within the 
2cr range. The relatively large deviations are probably due to the fact, that 
we use a non-regular triangulation, which may induce additional corrections 
to scaling. 

Assuming v = 1 we have next determined estimates for K c from the 
Binder-parameter crossings and the scaling of the various K max (L). The 
crossing points K x of the curves Ul(K) and UbL(K) with b > 1 approach K c 
as 

Kx =K c + K /(b 1 '"-l), (9) 

where k = k(L) does not depend on 6, and confluent corrections are ne- 
glected. This method, valid for large 6, is usually the most precise one and 
performing all possible least-square fits to (9) with L fixed we arrive at an 
estimate of the critical coupling K c = 1.023(1); see Fig. 1(a). Particularly in- 
teresting is that the crossings are close to their infinite volume value already 
on relatively small lattices. A similar observation with regards to Fisher ze- 
ros on spherical lattices was made in Ref . [9] . From the scaling of the various 
K ma , x , shown in Fig. 1(b), we can obtain further estimates of K c from linear 
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least-squares fits (assuming again v = 1). The error weighted estimate from 
all five sequences leads to K c = 1. 024(2), consistent with the crossing value. 
As our final value we used in further analyses 

K c = 1.023 ±0.001. (10) 

In particular we can now extract v also from the scaling of dU^ldK and 
d\n{\m\ p ) I dK at K c \ see Table 1. Here the errors reflect the combined un- 
certainties in K c and in the finite statistics. 

Also the values of Ul(K x ) vary very little with lattice size. Since in view 
of the statistical errors it was therefore impossible to perform an infinite- 
volume extrapolation (based on confluent corrections), we simply computed 
weighted averages over various subsets of the data. The results turned out to 
be quite insensitive to the precise averaging prescription, and as final value 
we quote 

U* = 0.55 ±0.01, (11) 

where the error is a very conservative estimate. Notice that this value is 
significantly different from the corresponding estimate for toroidal lattices 
(ss 0.610 — 0.615), i.e., it does depend on the global topology. 

To extract the critical exponent ratio 7/1/ we used the scaling x — 
L~'/"f 3 (x) at the previously discussed points of constant x, as well as the 
scaling of \' a ^ K c . The results are listed in Table 2. For the values of x 
at K c we could not find a fit of decent quality, however a fit at K = 1.020 
yielded 7/1/ = 1.756(6) with Q = 0.03. Similar difficulties with fits of \ at K c 
have been observed on regular lattices as well, see e.g. Ref. [10]. All values 
for 7/1/ are compatible with the Onsager value of 7/1/ = 1.75. The quality 
of the fit for Xmax can be inspected in Fig. 2. 

To extract the magnetical critical exponent ratio (3/u we used that {\m\) = 
L~P/ P fi{x) at all constant x-values. Another method is to look at the scaling 
of d{\m\)/dK = L^-^l v U(x). The fit results for (3/u and (1 - (3)/u are 
collected in Table 3. Using our average value for 1/u in Table 1 we obtain 
the final estimate of (3/u = 0.14(2). This result again is compatible with the 
Onsager value (3/u = 0.125. 

The specific-heat exponent a is numerically always the hardest quantity 
to estimate. For the case of the Onsager value a = we expect a logarithmic 
divergence like 

C(x,L) = a(x) + b(x)\nL. (12) 
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Indeed the data at the different fixed values of x can all be nicely fitted with 
this Ansatz. In particular, for the fit of C max with 10 data points shown in 
Fig. 3 we obtain a = 0.12(4), b = 0.361(9), with a total X 2 = 13.9 (Q = 0.09). 
We also did an unbiased three-parameter fit of the form 

C(x,L) = A(x) + B{x)L al ", (13) 

which gave us in the case of the fit of C max and 10 data points A = —5.6(15.), 
B = 5.8(15.), and a/v = 0.05(11), with a total X 2 = 13.8 (Q = 0.06). The 
barely improved \ 2 an( l l ne very small value of a/v, consistent with zero 
within error bars, clearly support logarithmic scaling; compare Fig. 3 where 
the two fits are indeed hardly distinguishable. The large errors on A and B 
can be understood by expanding L a l v = 1 + (a/v) In L + 0((a/ v) 2 ), showing 
that in the limit a/v — > the parameters of (12) and (13) are related by 
A + B = a and B = b/(a/v) — > oo. For comparison, we also included in 
Fig. 3 the best linear least-square fit with the KPZ prediction a/v = —2/3 
of the form 

C(x,L) = A(x) + B(x)L- 2/3 } (14) 

that resulted in A = 2.09(2), B = —6.9(2), and an unacceptable large total 
X 2 = 57.2 (Q = 1.6 X 10~ 9 ). Overall we can thus conclude that also for a/v 
our data is consistent with the flat-space Onsager value of zero. 

4 Concluding remarks 

We have performed a study of the Ising model on a spherical topology coupled 
to quantum gravity via the Regge calculus approach. Using the path integral 
measure JJdl/l we have found that, as on a toroidal topology, the critical 
exponents of the Ising transition agree with the Onsager exponents for regular 
static lattices, and the KPZ exponents are definitely excluded. The non- 
regular triangulation of the sphere seems to affect the finite-size behavior in 
a negative way, and one could probably obtain more accurate results, like 
those of Ref. [4] on the torus, by using a random triangulation of the sphere 
[11]. Unlike in the pure gravity case, where the global lattice topology enters 
in the formula for the string susceptibility exponent, it does not affect the 
critical exponents of the Ising phase transition. 
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Tables 



Fit-type 


AT 


1/u 


Q 


dU/dK at A^ ax 


7 


0.90(11) 


0.11 


dU/dK at K^ l) 


8 


1.10(6) 


0.44 


dU/dK at A^ ax 


8 


1.14(7) 


0.70 


dU/dK at At n f <|m|> 


8 


1.21(10) 


0.85 


at At" <m2> 


8 


1.23(10) 


0.87 


dU/dK at A" c 


8 


0.96(10) 


0.11 


d\n(\m\)/dK at A^ f <H> 


8 


1.13(5) 


0.53 


o?ln( m \)/dK at A" c 


8 


0.96(7) 


0.20 


dln(m 2 )/dK at A^ f <m2> 
dln(m 2 ) /dK at A" c 


8 


1.13(5) 


0.31 


8 


0.96(6) 


0.25 


weighted average 




1.08(5) 





Table 1: Fitting results for 1/u using K c = 1.023(1). Column N shows the 
number of fitted points, and Q denotes the standard goodness-of-fit param- 
eter. The average is computed by weighting each entry with its error. 



Fit-type 


N 




Q 


X at A^ ax 


10 


1.74(2) 


0.07 


X at KtT l) 
X at A^ ax 


10 
10 


1.739(8) 
1.745(6) 


0.66 
0.42 


X at At< |m|) 


10 


1.755(9) 


0.28 


X at Ktt 2) 
X' at K c 


8 

10 


1.73(2) 
1.79(3) 


0.29 
0.09 


weighted average 




1.744(6) 





Table 2: Fitting results for 7/1/ using K c = 1.023(1). Column N shows the 
number of fitted points. The average is computed by weighting each entry 
with its error. 
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Fit-type 


AT 


P/v 


Q 


(m) at A^ ax 


7 


0.07(4) 


0.17 


(m) at Kff 


8 


0.15(2) 


0.41 


(m) at A^ ax 


8 


0.15(2) 


0.34 


(m) at At n f <|m|> 


8 


0.20(5) 


0.87 


(m) at K^ m2) 


8 


0.20(6) 


0.71 


(m) at A" c 


10 


0.10(2) 


0.05 


average 




0.14(2) 








(l-P)/v 




d(\m\)/dK at A^ l} 
o?( m \}/dK at A" c 


8 
8 


0.93(3) 
0.89(4) 


0.27 
0.36 


average 




0.92(3) 





Table 3: Fitting results for (3/u and (1 — (3)/u using K c = 1.023(1). The 
number of fitted points is given under column N. The average is computed 
by weighting each entry with its error. 
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d=35 
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d=40 
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d=45 
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,1/v 



1/(b l,v - 1) (with v = 1) 



10 



K 



1.03 



1.02 



1.01 



1.00 



0.99 



0.98 



0.97 



b) 



-%f^ 



m . □ 



— e — 


max 




mf 


— 0- 


K z 

max 


- -X- 


- K ln(|m|) 

inf 


- - + - 


- T^-ln(m'm) 

inf 



0.02 



m — 



0.03 0.04 
L 1/v (with v = 1) 



0.05 



Figure 1: (a) The Binder-parameter crossings K x (b) with b = L'/L. As- 
suming v = 1, straight line fits through all data points yield K c = 1.023(1). 
(b) Finite-size scaling extrapolations of the pseudo-transition points K max vs 
L -1 /^, assuming v = 1. The error-weighted average of the extrapolations to 
infinite size yields K c = 1.024(2). 
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Figure 2: Finite-size scaling of the susceptibility maxima Xmax- The slope is 
compatible with the Onsager value 7/1/ = 1.75 for regular static lattices. 
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Figure 3: Finite-size scaling of the specific-heat maxima C max . Also shown 
are a logarithmic fit C max = a + 61nL, a power-law fit C max = A + BL a l v , 
and a constrained power-law fit assuming the KPZ prediction ajv = —2/3. 
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